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ABSTRACT 

A white dwarf (WD) gains substantial angular momentum during the accretion process that grows it toward 
a Chandrasekhar mass. It is therefore expected to be quickly rotating when it ignites as a Type la supernova. 
The thermal and shearing profile are important for subsequent flame propagation. We highlight processes that 
could affect the WD shear, during accretion as well as during the ^ 1000 years of pre-explosive simmering. 
Baroclinic instabilities and/or the shear growth of small magnetic fields provide sufficient torque to bring 
the WD very close to solid body rotation during accretion. The lack of significant shear makes it difficult 
to grow a WD substantially past the typical Chandrasekhar mass. Once carbon ignites, a convective region 
spreads from the WD's center This phase occurs regardless of progenitor scenario, and therefore it is of great 
interest for understanding how the WD interior is prepared before the explosive burning begins. We summarize 
some of the key properties of the convective region, which includes demonstrating that the mass enclosed by 
convection at any given time depends most sensitively on a single parameter that can be expressed as either 
the ratio of temperatures or densities at the top and bottom of the convection zone. At low Rossby numbers 
the redistribution of angular momentum by convection may result in significant shearing at the convective/non- 
convective boundary. 

Subject headings: accretion, accretion disks — convection — hydrodynamics — stars: magnetic fields — white 
dwarfs 



1. INTRODUCTION 

The fundamental role Type la supernovae (SNe la) play in 
determ ining the expansion history of the universe (iRiess et al.l 
|2004 and references therein) has added new urgency to un- 
derstanding the progenitors of these events. Theoretical work 
on modeling the flame dynamics of exploding white dwarfs 
(WDs) has demonstrated that the composition and energy of 
ejecta depend sensitively on the competition between flame 
propagation, instabi lities driven by turbulence, an d expansion 
of the WD (e.g., IHillebrandt & Niemeyed 12000). This has 
sparked new interest in understanding the evolution that oc- 
curs up until the onset of explosive burning. A critical aspect 
is the shear profile present within the WD. 

SNe la are expected to occur exclusively in accreting bi- 
nary systems. Therefore the role of angular momentum may 
be unique or at least especially important for these SNe. If ac- 
cretion drives shear in the WD, it may be an i mportant source 
of vis cous heating and material mixing. lYoon & Langed 
(l2004l) considered the accretion and evolution of a WD up 
until the point of ignition, including the effects of rotation. 
In their study angular momentum was transported primarily 
via the Kelvin-Helmholtz instabil ity, which lead to s ignif- 
icant shear throughout the WD. Saio & Nomotol (|2004) fo- 
cused on whether rotational effects can prevent accretion in- 
duced collapse (AIC) at high accretion rates (M = 3 x 10"^- 
10"^ Mq yr"'). They found that an AIC was not averted, but 
perhaps more interestingly, they also found the WD should 
be rotating nearly uniformly when the baro clinic instability 
is incl uded, a process that was neglected by lYoon & Langed 
(|2007). 

There is further opportunity for shear to develop during 
the carbon simmering phase. At carbon ignition, a convec- 
tive region grows at the WD center This envelops ^ 1 Mq 
over a timescale of 1000 yrs, until a burning wave com- 
mences. This phase has received increasing attention in recent 



years due to the realization that it sets t he initial temperatur e 
and density for the explosive burning (Lesaffre et al. 200^, 
as well as the distribution of ignition points (Woo slev et al] 
l200llWunsch & Wooslevl2004[lKuhlen et al.l2006h . Nuclear 
reactions during this time may change the neutron excess of 
the WD core, which is crucial for determining the ratio of ra- 
dioactive to non-radioactive nickel produced in the Type la 
explosion (Piro & Bildsten 2007b; Chamulak et al. 2001). 

In the present study we consider the opportunity for de- 
veloping shear during both of these stages. In ^ we re- 
visit estimates for the shear of accreting WDs and con- 
firm that the baroclinic instability limits the shear before the 
Kelvin-Helmholtz instability initiates. In addition we explore 
whether magnetohydrodynamic effects could reduce the shear 
even further. We conclude that the WD is nearly uniformly ro- 
tating at the onset of carbon ignition, and that viscous heating 
is negligible. In ^we explore the properties of convection 
present prior to explosive burning. The growth of this convec- 
tive region allows further opportunity for shearing. The three- 
dimensional nature of this problem makes it hard to defini- 
tively determine what occurs during this stage. We illustrate 
some general features expected for the interaction of spin and 
convection by appealing to observations and numerical exper- 
iments. In ^we conclude with a summary of our results and 
a discussion of future research. 

2. SHEAR PROFILE DURING ACCRETION 

Material accreted at a rate M reaches the WD surface with 
a nearly Keplerian spin frequency of JIk = (GM/R^y^^. The 
majority of the kinetic energy associated with this flow is dis- 
sipated in_j_boundary layer of thickness Hbl ^ R (as stud- 
ied bv iPiro & Bildstenll2004l) and never reaches far into the 
surface. Nevertheless, angular momentum is added at a rate 
MR^riji, so a torque of this magnitude should be communi- 
cated into the WD. It remains an open question whether all 
of this angular momentum is ultimately added to the WD. 
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If it was, the WD would reach overcritical rotation well 
befor e attaining the co nditions necessary for carbon igni- 
tion dMacDonaldl 1 1 9791) . Feedback with the accretion disk 
at high spin rates may provide a solution to this problem 
(jPaczynski 1991; Po pham & N arayan 1991). Observations of 
cataclysmic variables show a wide range of rotational veloci- 
ties, many of which are too low in comparison to the amount 
of an gular momentum that should have been accreted (iSioni 
[T999h . This may be related to the periodic mass loss th ese 
systems undergo in classical novae ( lLivio&Pringlelll998l) . 

For our present study the viscous mechanisms we consider 
have viscous timescales, fvisc = /v, where H = P/ pg is the 
pressure scale height, g = GMjr^ is the local gravitational ac- 
celeration, and V is the viscosity', that are much shorter than 
the timescale of accretion, face ^ 10^-10^ yrs. Due to this 
hierarchy of timescales we expect that the spin at any given 
moment can be broken into two contributions, which we write 
as 

l](r,f) = r!o(f) + Al}(r,r!o(f)). (1) 

The first piece, ^c(t), represents the solid-body rotation of the 
WD, which is increasing with time as accretion takes place. 
The second piece, Ail(r, ilo(0)i is the shear that must be 
present to transport the accretion torque into the WD. As long 
as the timescales obey f vise /face ^ 1, the viscosity is extremely 
efficient at transporting angular momentum. This means that 
the shear needed for transport is small, Afi <C ilo, and that 
the shear quickly comes into steady-state for a given 51()(f ), so 
that Ari does not explicitly depend on t. 

Since all of the transport mechanisms we examine oper- 
ate most efficiently in directions perpendicular to gravity (be- 
cause no work is performed), it is adequate to consider a 
structure composed of concentric spheres, each with con- 
stant 57. Transfer of angular m omentum is redu ced to a one- 
dimensional diffusion equation jFuiimotoll 19931) 
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We substitute equation ([T]i and assume that r is roughly inde 
pendent of time to obtain 
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The second term within the parenthesis on the left hand side is 
negligible since it is (9(Afi/51o) ^ 1. Multiplying both sides 
by Anpp- and integrating over radius we find 



I(r)——=4TTpvr -TT—, 
at or 



(4) 



where we have set the integration constant to zero to assure 
that the torque vanishes at the WD center and 



I(r): 



4Trpr dr 
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is the moment of inertia interior to a given radius. The rate 
of change of the WD spin is set by the accretion torque, so 
that Itotdflo/dt = MR-Qk, where /tot = I(R) is the total mo- 
ment of inertia. This gives us an equation for the shear 
a = dAn/dlnr, 



MR^n^Iir)/Iiat = ATipvr^a. 



(6) 



' This viscosity can represent either a molecular viscosity or a turbulent 
viscosity, but for the scenarios discussed in §2.1 a turbulent viscosity is ar- 
gued to be dominant. 



This resu h reduces to the ang ular momentum equation s pre- 
sented bv lFuiimotd ( Il993h and lPiro & BildstenI (l2007ah in the 
plane parallel limit by taking I{r) « /tot. With this expression 
we can now estimate the shear at a given depth as a function 
of ilo and M. The advantage of this approach is that VIq can be 
treated as a free parameter, and we are not restricted to follow 
the full history of angular momentum transport for the WD. 

Since for the remainder of ^ we assume that the approx- 
imations presented here hold true, we simplify notation by 
writing Vl as the solid-body rotation rate and a as the shear. 

2.1. Summary of Transport Mechanisms 

In the following sections we discuss some of hydrody- 
namic and magnetohydrodynamic instabilities expected to be 
present in the shearing WD core. We summarize three dis- 
tinct eddy diffusivities related to these instabilities and ex- 
plore their influence on the shear profile. This is not meant 
to be an exhaustive survey of all turbulent angular momentum 
transport mechanisms possible. Instead it is meant to high- 
light those that have been used in previous studies along with 
others that have received less attention in the past, but which 
we show to be important. 

2.1.1. Kelvin-Helmholtz Instability 

A hydrodynamic instability that has been popular for appli- 
cati qn to accretin g WDs is the Kelvin-Helmholtz instability 
(see lYoon & Langer,2004 . and associated work). This insta- 
bility is activated when Ri < 1 / 4, where the Richardson num- 
ber is 



N2 



Ri = 



and is the Brunt-Vaisala frequency, given by 



^ , dlnT\ 



(7) 



(8) 



(ignoring compositional gradients) where Q = 
-{dlnp/dlnT)p, Vad = (i91nr/c)lnP)ad is the adiabatic 
temperature gradient, and the asterisk refers to derivatives 
of the envelope's profile. The electron viscosity is the 
dominant molecular visco sity in the degenerate WD. S ubsti- 
tuting the viscosity from Nandkumar & Pet hickl (11984 ) into 
equation ^ results in a shear in which /?/ -C 1 / 4, and thus 
Kelvin- Helmholtz instability is expected. This confirms the 
result of lYoon & LangeJ (12004 ) that in the absence of other 
instabilities, Kelvin-Helmholtz instability is dominant, so 
that we view its associated shear as an upper limit. 

Kelvin-Helm holtz instability causes a turbulent eddy diffu- 
sivity given by (lFujimotolll993h 



(l-A-Ri)^/^ , 



(9) 



Since it only acts for Ri < 1 /4, it causes the shearing to evolve 
until /?/ = 1 /4 is marginally satisfied. Thus it is a good approx- 
imation to estimate the shear to be ctkh ~ 2A^ when secular 
instabilities due to thermal diffusion can be neglected (e.g., 
IZahnlfT992h . as is the case for the conductive WD core. Since 
electron degeneracy provides the dominant pressure, we ap- 
proximate the Brunt-Vaisala frequency as 



g ksT 
HZEp 



1/2 



(10) 
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where is Boltzmann's constant, Z is the charge per ion, and 
Ef is the Fermi energy for a degenerate, relativistic electron 
gas. We estimate for Kelvin-Helmholtz driven shear. 



(Tkh w 0.5 s gi 



1/2 .^1/2 

zl J73' (11) 
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where gm = g/lO'" cm s ^, /i^ is the mean molecular weight 
per electron, = p/lO'' g cm"-', and = T/IO^ K. This re- 
sult is in reaso nable agreement to the d etailed, time dependent 
simulations bv lYoon & Lange J (120041) . 

2.1.2. Baroclinic Instability 

Another hydrodynamic instability that may be important 
is the baroclinic instability (Fujimoto 1987, 1988; also see 
Gumming & Bildsten 2000). This instability arises because 
surfaces of constant pressure and density no longer coin- 
cide when hydrostatic balance is maintained under differen- 
tial rotation. In such a configuration, fluid perturbations along 
nearly horizontal directions are unstable, though with a suffi- 
cient radial component to allow mixing of angular momentum 
and material. When Ri is greater th an the critical barochnic 
Richardson number (lFujimotolll987h . 

GorioUs effects limit the horizontal scale of perturbations and 
the associated turbulen t visco sity is approximated from linear 
theory to be (iFuiimotdll993h 

1 Rise 



(13) 



3 /?/3/2 

In the WD interior Ribc 10, whereas the Richardson num- 
ber found for v = i/bc (using eqs. JS) and Q) is Ri ^ 10^, thus 
the limit Ri ^ Rise is clearly satisfied. 

The Richardson number can be estimated by combining 
equations (|6]l, d?), and (T3[ . giving 

= —^\r) q^nWY ^ ^ 

To illustrate how the shear depends on the properties of the 
interior we take 

3 M ^ fdynJ^K V/?/ ' 

where f^cc = M/M is the accretion timescale and fdyn = 
is the dynamical time. Using asc = N/Ri^^^ and /(r)//tot ~ 
(r/R)^, we estimate 



cbc ' 



10-4 I 



101 



0.1 



1/4 



(16) 



J V0-3S-1 

In ^2.2.11 we show that the timescale for viscous diffu- 
sion from the baroclinic instability is much shorter than the 
timescale over which accretion takes place. This means that 
this instability is able to quickly redistribute angular momen- 
tum as it accretes and the steady-state limit will be reached. 
Goupling this fact to the inequality (Tbc ~ o'kh/(2/?/) <C ctrh 
as demonstrated above, we conclude that the baroclinic insta- 
bility will limit the growth of shear long before the Kelvin- 
Helmholtz instabi hty can become active . This is consistent 
with the results of lSaio & Nomotd (l2004l) who included both 
the Kelvin-Helmholtz instability and the baroclinic instability 
and found nearly solid body rotation throughout the WD, even 
for accretion rates as high as 10"^ Mq yr"'. 



2.1.3. Tayler-Spruit Dynamo 

The last case we consider is magnetohydrodynamic insta- 
bilities, for which we apply the Tayler-Spruit dynamo (Sprui^ 
I2OO2I) . In this picture, shearing stretches any small compo- 
nent of poloidal magnetic field into a strong toroidal field. 
Once sufficiently large, the toroidal field initiates Tayler in- 
stabilities (non-a xisymmetric, pinch-lik e instabilities includ- 
ing stratification. lTayled 1 973l: ISpruit 1 999). which turbulently 
create poloidal field components that once again shear to be 
toroidal. This cycle continues and results in a steady-state 
field that transmits angular momentum via Maxwell stresses. 
In the limit when the magnetic diffusivity, rj, is much less than 
the thermal diffusivity, K, the minimum shear needed for this 
process to activate is (iSpruit,i2002i) . 



When a > (TTs.cnt, the effective viscosity due to the steady- 
state magnetic fields is 



1/4 



,3/4 



(18) 



We wait until the following section to present the shear asso- 
ciated with this process. 

2.2. Comparison Calculations 

We now calculate isothermal WD profiles and the shear pro- 
files for the mechanisms described above. Such models allow 
us to argue that solid body rotation is the most likely result. 

The WD models are computed by solving for hydrostatic 
balance with the effects of spin ignored since we generally 
consider spins of the order of 0. 1 Ok- We solve fo r p using the 
analytic equation of state from Paczyhski] (Il983l) . The impor- 
tance of Goulomb interactions is measured by the parameter 



1/3 



(ZeY Pa 
r= -— - = 35.7^-^ 



1/3 



(19) 



where A is the mass per ion and a is the ion separation. For 
the liquid phase, whe n 1 < F < 1 73, w e include the ionic free 
energy of Ghabrier & PotekhinI (Il998l) . The composition is 
set to 50% '^C, 48% '''O, and 2% ^^Ne by mass. 

In Figure [T| we compare the shear rates found by calculat- 
ing the effects of each viscous mechanism individually. In 
addition, we plot the critical shear required for the Tayler- 
Spruit dynamo to be activated, o-xs.crit (eq- iflTl ). All models 
use an isothermal temperature of 7} = 10*^ K and a WD mass 
of 1.37 Mq (i.e., near the critical mass for carbon ignition). 
We calculate (Tkh using the approximation ctkh ~ 2A^- Both 
(Tbc and ctts are calculated by substituting their associated 
viscosities into equation Q and solving for the shear using 
f7 = 0.mK = 0.67 s-i andM= 10-^Mq yr"'. 

Though the structure of an accreting WD is not exactly 
isothermal du e to the compe tition of compression and cooling 
(for example, Nomoto' 1982"), we can still make many conclu- 
sions using our simple model. Both ctkh and ctbc are consis- 
tent with what we estimated analytic ally. In fact, ctkh has a 
value within a factor of a few of what lYoon & Langed (|2004|) 
found for detailed accreting models (see their Fig. 7). Our cal- 
culations demonstrate that the baroclinic instability does not 
allow the shear to grow sufficiently for the Kelvin-Helmholtz 
instability to ever be important. Since ctbc ^ ^, the WD 
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Fig. 1 . — Comparison of the shear rate associated with three different vis- 
cous mechanisms as a function of mass coordinate: Kelvin-Helmholtz insta- 
bility (solid line), baroclinic instability (dotted line), and the Tayler-Spruit dy- 
namo (dashed line). We also plot the critical shear rate needed for the Tayler- 
Spruit dynamo to be activated (dot-dashed line). The WD is isothermal, with 
7; = 10* K,M= 1.37 MQ,Q = 0.1fiK = 0.67 s"', andM= 10"' Mq yr"'. 



should exhibit nearly solid body rotation. Also, note that this 
is for = 0.151k. Larger spins result in even smaller shear 
since <Tbc oc fl'^^^ (eq. |[T6l ). This is because at higher fl 
the relative amount of specific angular momentum in the disk 
versus the WD surface is smaller. 

We also consider the Tayler-Spruit dynamo in Figure[T] The 
thermal diffusivity is given hyK= 16f7sBrV(3 Cpup ), where 
CTSB is the Stefan-Boltzmann constant, Cj, is the specific heat 
capacity at constant pressure, and k i s the o pacity. For the 
opacity we include electron-scatteri ng jPaczvhsk i 1983), free- 
free, and conductive contributions (ISchatz et al.l ll999). The 
magnetic diffusivity is set a s 77 = TrkiTc ^ /(12e ^Kc), where 
Kc is the conductivi ty fro m ISchatz et al.l (Il999l) . The pre- 
scriptions given by Spruit (2002") imply an even smaller shear 
rate than the baroclinic instability. Associated with this small 
shear are steady-state radial and azimuthal magnetic fields, 
which we plot in Figure |2] This confirms that B,./B^ ^ 1, 
which is expected because the toroidal field growth is driven 
by shearing. The convection associated with the simmering 
phase (as studied in © may destroy this magnetic field for 
Mr < 1 Mq, but the field within the convectively stable re- 
gion will remain. Whether or not these fields are important 
for the later flame propagation when the WD is incinerated or 
for observations of the SNe la is an interesting question. Ac- 
cording to Piro_& Bildsten (2007a) both B,.,B<j) oc M'/- for the 
Tayler-Spruit dynamo. Thus any process or observational di- 
agnostic that is sensitive to the magnetic field strength would 
reveal something about M, an important discriminant between 
progenitor models. 

There is uncertainty in applying the Tayler-Spruit for- 
mulae to the case of an accreting WD. In the analysis 
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Fig. 2. — Estimates of (.solid line) and Br (dashed line) as implied 
by the Tayler-Spruit dynamo using the same WD model from Fig. [T] The 
fields associated with other accretion rates and spins can be estimated using 
oc M'/^O-l/s and B, oc A/'/^Q'/s iPh-o & Bildsten 2007a). Plotted for 
H/r < 1/2 is Bbuoy (dotted line), wh ich is the critical field needed for the 
buoyancy instability to occur (eq. 1211 ). 



presented by ISpruitI (l2002h . it is assumed that N > ^ > 
uja, where wa = B /[(4Trpy^^r] is the Alfven frequency. 
Such inequalities are appropriate for the radiative interior 
of the sun (for which this work was originally motivated). 
In the WD interior it is possible that > N since the 
B runt- Vaisala frequency is decreased by degenera cy effects 
(eq. IfTOl ). iDenissenkov & PinsonneauI3 (l2007h consider 
the effects of a large spin and find the effective viscos- 
ity of the dynamo is significantly reduced by a factor of 
(K/r^Ny^^(n/Ny^^(a/nf^^ <C 1. We hesitate from imple- 
menting their prescriptions because their results are based on 
purely heuristic arguments without the rigorous anal ysis of 
an ap propriate dispersion relation (as was provided in Spruit! 
12006). Since the baroclinic instability still contributes a large 
viscosity, our conclusion of solid body rotation is unchanged. 

An interesting possibility is that in the limit of large spin, a 
different instability other than Tayler instability is responsible 
for turbulently creating poloidal magnetic field components as 
is necessary for closing the dynam o loop. Th e magnetorota- 
tional i nstability ( Velikhov 1959; Chandrasekhar 1960; Fricke 
1969; AchesonH l978: Balbus & Hawley 1991, 1992; Balbus 
J995>) cannot provide closure to the dynamo since it requires 
dU,/dr < 0, which is opposite to what is found in the WD in- 
terior For similar reasons magnetic shear inabilities are also 
ruled out (Acheson 1978). Since the magnetic fields plotted 
in Figure |2] decrease with radius near the outer parts of the 
WD, it is possible that the m agnetic buoyancy instabi l ity oc- 
curs. Using the results from lAchesoiil (119781) . ISpruitI (119991) 
shows in the limit of cr/il -C 1 that such an instability arises 
for H/r< 1/2 when 



r] 



d\nB 
d\vir 



<0, 



(20) 



where we must remember that d\nB/d\vir < 0. This limit is 
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easiest to interpret if we solve for the magnetic field strength 
needed to satisfy this inequality, which is 



r/H-2y^^\d\nB/d\nr\ 



(21) 



This illustrates that the magnetic field must be sufficiently 
strong in comparison to the stratification (represented by A^) 
for the buoyancy instability to occur. We plot Z?buoy in Figure 
121 Since Bbuoy ^ B^, we conclude that the Tayler instability 
limits magnetic field growth and not the buoyancy instability. 
Apparently even when > N Tayler instability is the correct 
magnetohydrodynamic instability for closing the shear-driven 
dynamo in the core of accreting WDs. 

2.2.1. Viscous Timescales and Heating 

These WD models provide the local viscous timescale for 
angular momentum transport, fvisc = min[i/^,7?^]/z/, and the 
heating provided by viscous dissipation. In the top panel of 
Figure [3] we plot fyisc for the baroclinic instability and the 
Tayler-Spruit dynamo. Both are significantly smaller than 
the accretion timescale, thus our approximations presented at 
the beginning of ^ are justified. Compositional discontinu- 
ities present a possible barrier to angular momentum trans- 
port which we have ignored si nce we are focused on the WD 
core. iPiro & BildstenI (l2007ah show in the case of accreting 
neutron stars that such compositional changes can inhibit tur- 
bulent mixing, but do not alter angular momentum transport 
any more than introducing a slight spin discontinuity. Our 
assumption of steady-state transport in the core is therefore 
not affected. Figure [3] shows that there is a clear hierarchy 
of timescales. If both mechanisms were acting, the Tayler- 
Spruit dynamo acts sufficiently rapid that it limits the shear 
before the baroclinic instability becomes important. 

The viscous heating per unit mass is 



Evisc = Tl^CT 



SO that the total luminosity is 



(22) 



(23) 



In the bottom panel of Figure |3] we plot the integrand of this 
expression, dLyi^c/dMr = eyisc, for the baroclinic and Tayler- 
Spruit cases in units of L© Mq . We denote the heating in 
this manner to emphasize that with these units one can easily 
"integrate by eye" to determine the total luminosity over any 
region of the WD. The process of accretion and compression 
leads to additional heat ing, which has been studied analyti- 
cally in Appendix B of iTownslev & BildstenI ( 120041) . Using 
their results, we approximate that within the degenerate core. 



dL. 



'comp 



3 AtbT M 
dMr 5 ^tnip Mr ' 



(24) 



where /x, is the mean molecular weight per ion. This is also 
plotted in the bottom panel of Figure [3] {dashed line), which 
shows that the viscous heating is dwarfed by this compres- 
sional heating. Not plotted on Figure [3] is the compressional 
heating in the non-degenerate envelope. This total integrated 
contribution to the heating is larger than t hat of the de generate 
core for the accretion rates of interest (Nomotol ll982h . 

We note that Figure [3] shows a model with M = 
10"^ Mq yr"', and that at shallow depths baroclinic heating 
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Fig. 3. — In the top panel we compare the timescales for viscous trans- 
port of angular momentum via the baroclinic instability (solid line) and the 
Tayler-Spruit dynamo (dotted line). The bottom panel shows the rate of en- 
ergy generation for both viscous mechanisms, as well as t he heating expected 
from the compression of the WD from accretion (eq. 1241 , dashed line). The 
WD models are the same as described in Fig. [T] 



is merely a factor of a few le ss tha n compressional heating. 
Using the scalings derived in ^2.1.2l we find that dL^c / dMr 
M^/"*, whereas dL^omp/dMr cx M. Thus viscous heating from 
baroclinic instabilities is important when M > 10"^ Mq yr"', 
which was demonstrated in some of the models considered by 
iS_aio&Nomoto(2004). 

3. THE CONVECTION PHASE AND REDISTRIBUTION 
OF ANGULAR MOMENTUM 

The explosive burning that incinerates the WD in the SN 
la does not begin right at carbon ignition. Instead, the 
carbon simmers f or ^ 1000 yrs , which we summarize in 
Figure |4] (also see 'Woosl ev et ani2004t IWunsch & Woosleyl 
2004; Lesaffre et al. 200g). As this stage is important 
for setting the initial conditions for explosive burning, 
there has recently been significant progre ss in modeling 
it numerically ( Hoflich & Stein 2002; Ku hlen et alj l2006t 
IStein & Whederi i2006t lAlmgren et_aL| l2006allbl) . For our 
purposes here we present simpler models that nonethe- 
less capture many of the crucial features of this stage. 
For simplicity we igno re th e convective Urea process 
(Paczvhski 1972; Bruenn 1 19731; ICo uch & Arnett 1975; Ibe^ 
1978a b, 1982; Barkat & Wheeler 'l 990; Mochkovitch 199^ 
Stein e t al. 1999; Bisnovatyi-Kogan 2001; Lesaffre et al] 
I2OO5I: IStein & Wheeler) 12006) . Ignition first occurs when 
heating from carbon burning beats neutrino cooling. This 
is plotted as thin solid lines that show the range of theoret- 
ical ignition curves expected for a uniform composition of 
'^C and '^O, with th e middle line giving the optimal result 
(lYakovlev et alJl2 006). 

Simmering is followed by calculating a series of hydro- 
static WD models each with a different central temperature. 
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mal with Ti= 10** K. The dotted lines show characteristic heating timescales 
as labeled. 



is the convective overturn timescale, Vconv is the characteristic 
convective velocity, and we have used a mixing length equal 
to the minimum of either the pressure scale height, H, or the 
radius of the convective core, /?conv The exact condition re- 
quired for the deflagration to begin depends on what fraction 
of the convective zone is still responding to the energy input 
and on th e extreme temperatu re sensitivity of c arbon burning 
e PC r ^^, [Wooslev et al.ll2004h . For this reason lUesaffre et all 
20061) explore f^^n,, = atu. where a < 1 parameterizes this un- 
certainty. Since we are only roughly concerned with resolving 
the end of the simmering phase, we take fconv ~ h to find that 
simmering should end when ?/, w 7 s at a central temperature 
and density of 7; k. 7.8 x 10^ K and p, « 2.6 x 10^ g cm'^ 
This is roughly in ag reement with the estimates presented by 
'Wooslev et al. ( 2004) using the Kepler stellar evolution code 
(^ea_ver et al. 1978). 

The series of events that take place during the simmering 
phase are expected to occur regardless of the progenitor sce- 
nario under consideration as long as ignition occurs near the 
WD center. Therefore, it is important to study any processes 
that may influence the subsequent explosive burning of the SN 
la. In the following sections we describe some of the main 
characteristics of the convective region, with a focus on its 
mass and the velocities of convective eddies. Afterwards we 
estimate how angular momentum may be redistributed by the 
convection. 



but with a fixed mass of 1.37 Mq. We use the same micro- 
physics and composition as in 32.21 and ignore the compo- 
sitional changes due to the burning. The thick dashed line 
on the right of Figure H] traces out the central temperature 
as a function of density as the simmering takes place. The 
density decreases as the WD is heated and expands. The 
thick solid lines show example profiles demonstrating how 
the convection grows as the central temperature increases. 
Since the thermal conduction timescale is so long (~ 10^ yrs), 
the convection is very efficient and the convective profile is 
nearly an adiabat. The adiabatic index depends sensitively on 
the fact that the ions are in a highly correlated liquid state 
dTownslev & Bildstenll2004l) . 

/^^lnr\ 0.91+0.14ri/3 

\d^p) 1.22 + 0.4iri/3' ^ ^ 

where F is the parameter defined in 32.21 For the range of 
models we consider, we find (d\aT ld\vi p)^^ « 0.49-0.52 
(roughly consistent with I Wooslev et al.ll2004l) instead of the 
usual Ki2/3> that would be appropriate if the ions were an 
ideal gas. The top of the convective region is set by when 
it reaches the isothermal temperature of 7} = 10'^ K, in other 
words, whe re the it's entropy matc hes that of the initial ther- 
mal profile (IHoflich & Steinll2002h . 

Also plotted in Figure ID are dotted lines that show char- 
acteristic values for the heating timescale, f/, = CpT jt, where 
e is the e nergy generation rate for carbon burning from 
ICaughlan^& Fowler ( 198 81) with s trong screening included ac- 
cording to [Salpeter & van Horiil (IT969I) . This gives a rough 
estimate of the time until the burning wave commences. For 
the majority of the simmering stage, the true timescale to heat 
the convective region is significantly longer because the heat 
capacity of the entire convecti ve zone much be ta ken into ac- 
count (see the discussion in .Piro & Chang||2008l) . The defla- 
gration begins once th ~ fconv, where 

fconv = min[//,/?(;onv] /Vconv, (26) 



3.1. Properties of the Simmering Core 
3.1.1. Convective Core Mass 

The boundary between the convective and non-convective 
zones reaches out to densities as shallow as ^ 10''- 10^ g cm"^ 
(Fig. m and encompasses a mass of Mconv ^ 0.9- 1.2 Mq. 
Even though the exact transitional depth depends sensitively 
on the WD's isothermal temperature as wel l as the com- 
positi on and microphysics near the boundary (iPiro & ChangI 
l2008h . it is useful to estimate Mconv 

In order to gain some intuition for the expected dependen- 
cies we solved for Mconv using the Lane-Emden equation for 
a poly trope of index 3 in the limit of small radius. This 
showed us that Mconv depends on only two parameters: the 
mean molecular weight, pe, and the ratio of the isothermal 
and central temperatures, Ti/Tc- 

Unfortunately, this analytic result cannot be directly applied 
to the problem at hand to estimate Mconv (because the small 
radius limit is not accurate). We therefore search for the true 
dependence numerically. In Figure |5] we plot Mconv versus 
Ti/Tc for four different models that sample a range of central 
densities and isothermal temperatures. This shows that the 
relation is indeed very tight. A fit to the solid fine in Figure |5] 
gives 

Mconv = 1.48 Mq (^^^ (^1-1.20^^. (27) 

This result could also be written as a ratio of densities since 
the convective zone follows an adiabat with approximately 
p oc T^. Changes in p^ are primarily due to ^^Ne, which 
gives a range p^ ~ 2.00-2.01 for a mass fractions range 
X(^^Ne) = 0.00-0.06. Therefore pe can change Mconv by less 
than 1%. We note that an important complication we have 
omit ted is the presence of gradients in the C/O abundances 
(e.g. Lesaffre et al. 2006). 

This result (eq. ll27l ) is maybe not too surprising. It is well 
known that for a relativistic, degenerate equation of state that 
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Fig. 5. — The mass of the convective core, Mjonv, as a function of the ratio 
of the isothermal and central temperatures, Tj/Tc. We consider models with 
a range of central densities, pc, and isothermal temperatures, Tj, as labeled 
within the plot. All models trace out a similar relation, which demonstrates 
that Mconv = Mcom(Ti/Tc) is a good approximation. 

the characteristic mass of a self-gravitating object (the Chan- 
drasekhar mass) is independent of the central density and only 
depends on /i^- In addition, nuclear reactions during the sim- 
mering may also complicate things, which we have ignored. 
Nevertheless, it is an important relation to keep in mind. It 
tells us that independent of the complicated previous history 
of a SN la progenitor during the simmering phase, at any 
given time Mconv depends most strongly on just a single di- 
mensionless number. 

3.1.2. Convective Velocities 

Besides the size of the convective zone, another important 
property is the speed of the convective eddies. Using esti- 
mates from mixing-length theory, the characteristic convec- 
tive velocity, Vco nv, is related to Fmnv for efficient convection 



via dHansen & Kawaleri[T994l) 



CnT n 



QgH' 



= pvL 



(28) 



where we have used a mixing-length that is equal to the scale 
height, H. If the convective region was static during simmer- 
ing, then Fconv would simply be equal to the integrated nuclear 
energy generation. Instead the convection zone is growing 
and energy must be expended to heat new material. Also, the 
convective flux must be nearly zero at the convective bound- 
ary because of the long conductive timescale here 10* yrs). 
These effects cause a decreased /^onv in the outer parts of the 
convective regi on, which we inves tigate in more detail in a 
separate study ( Piro & Changll2068h . We borrow these results 
for presenting Vconv here. 

In the upper panel of Figure |6] we plot Vconv as a function 
of mass coordinate using equation (|28] l. The characteristic 
value of Vconv ~ 10^ cm s~' during the late sta ges of convec- 
tion is consistent with the analytic estimates of IWooslev et al.l 




0.001 



Fig. 6. — In the top panel we plot the characteristic convective velocity, 
Vconv , as a function of mass coordinate for various times during the simmeiing 
phase. Each line is labeled by the heating timescale, f;, = CpT/e, which is 
meant as a rough estimate for the time until the deflagration begins. In the 
botto m p anel we plot the convective eddy turnover frequency, ficonv = 'coiw 
(eq. 1261 ). for the same four models. The dashed lines denote a reasonable 
range for the WD spin U (0.1 - 1.0)r2K. 



( |2004|) . The turbulent conditions prepared in the core will 
be important for understanding the subsequent propagation of 
bubbles and burning. For example, Zingale & Dursi (200'^ 
show that the turbulence may act as a viscous drag with the 
largest effect on the smallest bubbles. Another interesting fea- 
ture is that the material above the convection is devoid of this 
turbulence. The burning properties could change significantly 
as a flame passes into this relatively "quiet" region. 

The importance of the WD spin for the convection is deter- 
mined by the convective Rossby number 



(29) 



where fir 



(see eq. 



) is the eddy turnover fre- 



quency. In the bottom panel of Figure |6] we plot ilconv along 
with the expected range of spins relevant for SN la progenitors 
(dotted lines). For the entirety of the simmering stage we find 
Ro <l, showing that spin should have a non -negligible influ - 
ence on the convection (for example, see Ku hlenetal.ll200^ . 
This will be an important consideration in the following sec- 
tions. 

3.2. Rotation Profiles during Simmering 

We next consider how the spin changes as convection grows 
in the core. For the 10^-10^ yrs that accretion takes place 
we have shown that hydrodynamic or magnetohydrodynamic 
instabilities severely limit shearing. During the last stages 
of simmering the WD is changing on the timescale of hours 
to minutes, so these viscosities will have less influence and 
stronger shearing is possible. 

Given the conclusions made in ^ we consider uniform 
rotation to be a reasonable starting point for when carbon 
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first ignites. Heating and expansion generally causes the WD 
to spin down (like the freshman physics figure skater prob- 
lem). In the following calculations we assume that angular 
momentum is conserved locally on spherical mass shells in 
non-convective regions. The mixing by convective eddies is 
important for determining the spin within the convective core. 
Such effects are difficult to estimate, so we consider two pos- 
sible rotation laws. 

3.2. 1 . Uniform Spin within the Convection 

The first rotation law we explore is that convection enforces 
solid body rotation. This is the simplest and mos t nai ve as- 
sumption one can use. The study by Kumar et"al] (1199 5) ar- 
gues that such a rotation law is reasonable when convective 
eddies scatter elastically. Unfortunately, it is not clear whether 
such a picture of the eddies as separate fluid elements is ap- 
propriate. Elastic scattering is in fact counter to the standard 
assumption in mixing-length theory that eddies deposit their 
entire angular momentum and entropy into the ambient fluid 
with each scattering. This result also assumes an isotropic 
scattering whereas rotation clear l y brea ks this symmetry. In 
the anisotropic case iKumar et alJ (Il995h find that in principle 
angular momentum can be transferred outward by convection, 
an idea we consider in more detail in the next section. 

In Figure |7] we plot the evolution of the WD spin during 
the simmering stage assuming solid body rotation within the 
convective zone. The initial spin is taken to be uniform with 
ft = 0.67 s"' (dashed line). The overall trend is spin-down 
from heating and expansion. Uniform rotation in the convec- 
tive zone leads to shear at its top with a discontinuous velocity 
of Ay 10^-10*' cm s"' at late times. 

3.2.2. Outward Transport of Angular Momentum by Convection 
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Fig. 8. — The spin along the equator versus mass coor dinate assuming 
a rotation law within the convective zone give be eq. j30t , using the same 
background models as Fig. [7] 

The other rotation law that we consider is motivated by nu- 
merical simulations and observations of the Sun. In theoreti- 
cal and numerical studies where /?o < 1, a generic feature is 
that angular momentum is transported outward away from the 
poles and t owards the equ ator ( Oilman 1979; Miesch 200C)I; 
iBrun & fbo mre 2002; Br owning et alJl2004l) . This appears to 
be controlled by t he largest scale eddies that are most influ- 
enced by rotation dBrowning et al.ll2004l) . Such features are 
qualitatively consistent with helioseismic measuremen ts that 
map t he outer convective region of the Sun (Thompson et alj 
I2003L and refer ences therein). The main disparity is that the- 
ory and numerics generally result in a more Taylor-Proudman 
like spin profile whereas the Sun's velocity contou rs are more 
radial . This may be due to the tachocline (Miesc h et alj 
|2006|) . which is present in cases with surface c onvection. The 
study o f core convection in rotating A-stars bv lBrowning et alj 
(120041) . which also finds a Taylor-Proudman spin profile, may 
be the most relevant comparison to the WD case. 

To mimic the general features of the spin profiles described 
above we consider the following rotation law. 



rsinc 

^rnnv 



+ 1 



(30) 



where 0^ is the spin at the WD's center, is the fractional 
change in spin across the convective zone, and d is the lat- 
itude measured from the pole. For calculational purposes, 
we choose /5 = 0.4. This is in reasonable agreement with 
the results o f Browning et al. (2004) and solar observations 
([Thompson et al. 2003), which both have a similar value of 
Ro ^^.\ in comparison to our case here. 

In Figure [S] we plot the spin profile along the equator using 
equation ( l30l ) as the rotation law within the convection. The 
outward transport of angular momentum is much more promi- 
nent in comparison to the uniform rotation case considered 
earlier Since the rotation is cylindrical, at higher latitudes the 
shear is significantly smaller. The typical velocity jump at the 
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Fig. 9. — The heating available from the sheai' at the convective/non- 
convective boundary as a function of the boundary's mass coordinate. These 
are found for both of the rotation laws we consider for the convection of solid- 
body rotation (solid line) and outward angular momentum transport (dotted 
line). As a comparison we show the thermal energy at the top of the convec- 
tion zone, CpT (dashed line). 

equator is AV 10''- 10^ cm s"', which may be comparable 
to the speeds of the burning fluid elements that will be buoy- 
antly rising through the WD once explosive carbon ignition 
occurs. This could in principle shear out the burning, enhanc- 
ing it because of the increased surface area. 

3.2.3. Heating at the Convective Boundary 

The shear present at the convective/non-convective bound- 
ary represents free energy that has been made available be- 
cause of angular momentum transport by convection. If this 
shear persists throughout the convective stage, then it may 
have interesting effects on the subsequent flame propagation. 
The shear is not smeared out by convective overshooting. Us- 
ing an overshooting length of L « V^^^^/ilg) and a sound 
speed Cj, we estimate that L/H ^ (Vconv/cj)^ <C 1, since the 
convective velocities are very sub-sonic. On the other hand, 
viscous processes may act at this interface, damping out the 
shear and leading to heating. Assuming all of the shear is con- 
verted to heat gives a local heating per unit mass of ^ (AV)^. 
This is basically an upper limit of the energy generation rate 
at any given time. 

In Figure |9] we plot the free energy per unit mass available 
from the shear at the convective/non-convective boundary for 
both rotation laws. This is shown to be less than the thermal 
energy at the same location, CpT. We note that these results 
are for an initial spin of « 0.1 fix (chosen so that the spin 
would have minimal impact on hydrostatic balance). Since 
Ay oc fl, the heating can be up to « 100 times larger than 
what we estimate here. In the case where angular momen- 
tum is transported outward by convection, dissipation of shear 
may be able to heat the outer layers of the star This may be an 
important effect for priming the thermal profile and expanding 
the WD before the flame propagates past. 

4. DISCUSSION AND CONCLUSION 

In this study we have investigated the internal rotation pro- 
file for SN la progenitors, focusing on the stages of accretion 



and simmering. Since these progenitors are formed exclu- 
sively in binary systems, they contain a enormous amount of 
angular momentum. This is important for the WD structure 
and the subsequent flame propagation during the explosive 
burning of a SN la. 

During the accretion phase we compared three different 
mechanisms for angular momentum transport. Either the 
baroclinic instability or the Tayler-Spruit dynamo limit the 
shear to values so low that the Kelvin-Helmholtz instability 
cannot occur. We argue that the WD will have nearl y uniform 
rotation, consistent with the results oflS aio & Nomo to (2004). 

This conclusion has important implications for ex- 
tremely energetic S Ne la such as SN2003fg (SNLS-03D3bb, 
iHowell et alj l2006h . which has been argued to be from a 
super-Chandrasekhar WD dJeffrev et alj 1200 6). One well- 
known way to get a mass this high is v ia strong differ- 
ential rotation dShapiro & Teukolskvl I1983L and references 
therein). If uniform rotation is required then the WD mass 
can only be increased by a few percent over the normal Chan- 
drasekhar limit. We should also mention that just because 
the WD mass is higher, the chemical yields of the deflagra- 
tion or detonation of a strongly differentially rotating WD 
will not ne cessarily match the high nickel mass needed for 
SN2003fg dSteinmetz et al . 1992; Pfannes 2006). For this rea- 
son, other interpretations have been suggested for this SN la 
(Hillebrandt et al. 2007). 

The simmering phase is another opportunity for differential 
rotation to develop within the WD. This occurs due to heating 
and expansion as well as angular momentum mixing by con- 
vection. The full impact of the convection is hard to determine 
with certainty. Given the characteristic convective velocities, 
we expect Ro < I, thus the influence of rotation cannot be 
discounted. We consider two rotation laws for the convective 
zone, a uniform rotation case and another case where the spin 
increases outward. In either case the most dramatic shear- 
ing appears at the convective/non-convective boundary. This 
shearing may have an impact on the propagation of flames or 
bubbles, and in some extreme circumstances may hold enough 
free energy to alter the thermal profile of the non-convective 
envelope. 

Given the impact rotation could potentially have for SNe la, 
the time appears ripe for including such effects in simmering 
and flame pro pagation calcula t ions. S imulations such as those 
performed by lBrowning et al.l (l2004h are not yet able to reach 
the characteristic fluid parameters expected during the sim- 
mering stage (such as the Reynolds and Rayleigh numbers; 
see th e summary in lWoosley et al.ll2004l : IWunsch & Woosleyl 
l2004l) . Nevertheless, « 0.1 is well within reach, as has 
been shown by the simulations of iKuhlen et alj d2006h . The 
biggest difficulty for the future will be following the secular 
evolution of the WD as the convective zone grows. Since the 
carbon burning is so temperature sensitive, the majority of the 
evolution takes place at late times (see Fig. |2]and[8]l. This will 
help limit the duration that must be simulated. 
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